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Abstract 



The Galactic orbital parameters of 159 cataclysmic variables in the Solar neigh- 
bourhood are calculated, for the first time, to determine their population types using 
published kinematical parameters. Population analysis shows that about 6 per cent 
of cataclysmic variables in the sample are members of the thick disc component 
of the Galaxy. This value is consistent with the fraction obtained from star count 
analysis. The rest of the systems in the sample are found to be in the thin disc 
component of the Galaxy. Our analysis revealed no halo CVs in the Solar vicinity. 
About 60 per cent of the thick disc CVs have orbital periods below the orbital period 
gap. This result is roughly consistent with the predictions of population synthesis 
models developed for cataclysmic variables. A kinematical age of 13 Gyr is obtained 
using total space velocity dispersion of the most probable thick disc CVs which is 
consistent with the age of thick disc component of the Galaxy. 

Key words: 97.80. Gm Cataclysmic binaries, 98.10.+Z Stellar dynamics and 
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1 Introduction 



Cataclysmic variables (CVs) are short-period interacting binary stars, consist- 
ing of a white dwarf primary and a low-mass late spectral type secondary star. 
The secondary star fills its Roche lobe and transfers matter to the primary via 
a gas stream and an accretion disc. Accretion disc formation is prevented in 
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systems with strongly magnetic white dwarfs in which mass accretion contin- 
ues through acc retion columns . For a detailed review of observational proper- 
ties of CVs, see IWarnerl (119951 ). 



Stehle et al.l (119971 ) studied the long-term evolution of cataclysmic variables 
as a function of the secondary star metallicity. They showed that Pop II 
CVs with a low metallicity secondary star have a detac hed phase with a 



smaller orbital period width, a shorter minimum period (jPaczynskil . Il967l ) 
and a slightly higher mass transfer rate, resulting in shorter evolutionary 
timescales compared to CVs where the secondary star has a Solar chemi- 
cal composition. According to their population synt hesis model, most Pop I 



CVs are expected to be found below the period gap (IVerbunt fc Zwaanl. 11 981 



Rappaport et al.| . |l982 . 1983 ; Paczynski fc Sienkiewicz . 1983 ; Spruit fc Ritter 



19831 : lKingl . ll988[ IShao fc Li I . l2012h . IStehle et all (Il997l ) express that the high 
7 velociti es (systemic velocit i es or center of mass velocities) of some systems 
found by Ivan Paradijs et al.l (119961 ) suggest these systems to be Pop II CVs. 
However, most of these CVs are magnetic systems (DQ Her and AM Her 
stars), where the Doppler-shifts of spectral lines originate mainly from the ac- 
cretion stream. Consequently, the errors in the 7 velocities may be noticeably 
high. Still, finding magnetic systems below the period gap should not be a sur- 
prise as these systems con centrate towar ds shorter orbital periods, with little 
evidence for a period gap (jWarnerl . Il995l ). Interestingly, this narr ow (or none) 
period gap is consistent with th e predicti o ns fro m the study of IStehle et al. 
( 119971 ) . It should be noted that lAk et al.l ( 120101 ) found magnetic systems to 
be much older than non-magnetic systems while they also emphasized doubts 
about the reliability of 7 velocities obtained from the observations of magnetic 
systems. 



the Galaxy ( 


Hawkins & Veron, 


1987: 


Howell. Mitchell & Warnock III. 


1987; 


Howell & Szkodvl. 


1990; 


Drissen et al.. 


1994; 


Sheets et al.. 


2007; 


Uthas et al.. 


2011; 


Imamura & Tanabd. 


20121). Clearly, it is hard to detect thick disc and 



halo CVs if they are not outbursting systems. CVs in globular clusters could 
be detected from their outbursts, emission lines, X-rays or their blue colours. 
Low metallicity values ex tracted from spect ra can be another indicator for 



thick disc and halo CVs ( IStehle et al.l . 119971 ) . However, spectroscopic confir- 
mation of these detections is difficult as they are very faint systems. For a 
detailed discuss ion on the detection of CVs in globular clusters, we refer to 



Knigge's ( 1201 ll ) review and references therein. 



A reasonable observable sample of Pop II CVs (thick disc and halo CVs) 
can only be f ound at vertical distances from the Galactic plane z > 2 kpc 



( IStehle et al.l . 119971 ) . However, the faintness of CVs restricts most of their 



photometric and spectroscopic studies to the Solar neighbou rhood. Although 



some abundance anomalies in UV and IR were found for CVs ({Hamilton et al. 
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201ll ; lGansicke et all 120051 ). there are no reliable metallicity measurements for 
CVs. In the absence of metallicity measurements, only kinematical and the dy- 
namical methods are expected to provide reliable results in the recognition of 
the thick disc and halo CVs in the Solar neighbourhood. Although the de- 
velopment of sensitive CCD cameras, spectrographs and larger size telescopes 
made it possible to observe relatively faint CVs, we lack the necessary obser- 
vational information to understand the kinematical properties of CVs in the 
Galactic scales. 

In the next section we use the kinematical properties of 159 CVs in the solar 
neighbourhood to distinguish different Galactic populations among them. 



2 The data 



2.1 Input Parameters 



The kinematical data used in this study are taken from lAk et al.l (120101 ) . The 
most important inputs used in their study are distances, 7 velocities and 
proper motions. The proper motions of C Vs were mostly obtained from the 



NOMAD Catalogue (IZacharias et all 120051 ). The ty pes, equatorial c oordi nates 



and orbital p eriods of CVs were m ostly taken from 
tion 7.7) and lDownes et al.l ( 1200 ll ). 



Ritter fc Kolbl (120031 . Edi- 



The distances were predicted using the PLCs relation of lAk et al.l ( 120071 ) . The 
PLCs rel ation is based on orb ital period and Two Micr on All Sky Surve y 



(2MASS, ISkrutskie et al.l . |2006|) JHK S photometric data flCutri et al.l . |2003|) 



This relation is reliable and valid in the ranges 0.032 <P(d) < 0.454, —0.08 < 
(J ~ H) < 1.54, -0.03 < (H - K s ) < 0.56 and 2.0 < Mj < 11.7 mag. 
For a d etailed descr iption of the method by the PLCs relation, we refer to Ak 
et al.'s (120071 120081 ) studies . The distances obtained from this relation differ 
i n general less than 4% fro m those obtained from trigonometric parallaxes 
(iGarietv fc Ringwaldl . lioiih . 



The last important inputs used in lAk et al.l (120101 ) are the radial velocities 
which are used to calculate total space velocities. The radial velocity with 
respect to the Sun comes from measurements of Doppler-shifts in spectral 
lines. However, one has to take here into account that CVs are binaries, and 
thus on e has to determ ine the radial velocity of the cen ter of mass of the 
system . lAk et al.l (120101 ) adopted the criteria defined by Ivan Paradijs et al. 
(119961 ) when collecting 7 velocit ies from the literature an d merged these 7 
velocities with those collected by Ivan Paradijs et al.l (119961 ). Radial velocities, 
and consequently 7 velocities, derived from emission lines are likely affected 
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by the motion in the accr etion disc or the matter stream falling on the disc 
from the secondary. Thus, lAk et al.l (120101 ) analyzed 7 velocities statistically 
and looked for possible systematic errors in the 7 values obtained from emis- 
sion lines. They concluded that there is no substantial systematic difference 
between systemic velocities derived from emission and absorption lines and 
that the o bserv ed 7 velocities can be reliably used for statistical analysis (see 
Ak et all , bold , for details). 



From the celestial coordinates (a, 5), prop er motion compo nents (fi a cos S, fis), 
systemic velocity (7) and the parallax (7r). lAk et al.l (120101 ) computed Galactic 
space velocities and their propagated error s with respect to the Sun usin g the 
algorithms and transformation matrices of I Johnson Ik Soderbloml (119871) . Al- 
though the sampled CVs are relatively nearby objects, Ak et al. ( 2010[ ) applied 
corrections for differentia l Galactic rotatio n to space velocities as described in 
Mihalas fc Binnev! fjl98lh . lAk et all (l2010f ) analyzed the propagated errors of 
space velocities, with respect to LSR (Local Standart of Rest), and refined 
their sample by removing systems with a total space velocity error S err > 30 
km s _1 . Although they collected input data from the literature for 194 CVs, 
this analyses decreased the number of usa ble system s in th eir sample to 159. In 
this study, the final sample of 159 CVs in lAk et al.l (120 10f ) is used to calculate 
the Galactic orbital parameters of systems. 



2.2 Calculation of Galactic Orbits 



In order to determine possible Galactic orbits of CVs, we first perform test- 
particle integration in a Milky Way potential which consists of a logarithmic 
halo of the form 



$ h aio(r)=u 2 m(l + -J, (1) 

with vq = 186 km s~ x and d — 12 kpc. The disc is represented by a Miyamoto- 
Nagai potential: 



<IW(/?-:) (2) 

& + (a d + J? + b\ 



with Md = 10 11 Mq, ad = 6.5 kpc and = 0.26 kpc. Finally, the bulge is 
modelled clS cL Hernquist potential, 



$bui g c(r) = -^p^, (3) 
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using Mb = 3.4 x 10 10 M Q and c = 0.7 kpc. The superposition of these 
components gives a good representation of the Milky Way. The circular speed 
at the Solar radius is taken ~220 km s -1 . The orbital period of the LSR is 
Plsr — 2.18 x 10 s years while V c = 222.5 km s _1 denotes the circular rotational 
velocity at the Solar Galactocentric distance, Rq = 8 kpc. The same formulae 
were alread y used to determine the Gala ctic orbits of obj ects from different 



classes, e.g. ICoskunoglu et al.l (120121 ) and iBilir et al.l ( 120121 ) 



In order to analyse Galactic orbits of CVs, the mean radial Galactocentric 
distance (R m ) is taken into account as a function of the stellar population and 
the orbital shape. We consider the planar and vertical orbital eccentricities, 
e p and e v , respectively. R m is defined as the arithmetic mean of the final 
perigalactic (R p ) and apogalactic (R a ) distances, and z max and z min are the 
final maximum and minimum distances, respectively, to the Galactic plane 
whereas e p and e v are defined as follows: 



e - ^ 5^ (a) 

t^a + -ftp 



and 



\Zmax \ ~\~ pmtn 



R,, 



(5) 



respectively, where R m = (R a + R p )/2 (jVidojevic fc Ninkovid . 120091 ). Due to 
^-excursions R p and R a can vary, however this variation is not more than 
5%. Calculated orbital parameters of the 159 CVs are listed in Table 1. The 
columns of the table are organized as follows: name, equatorial (a, 5) coordi- 
nates, the final perigalactic (R p ) and apogalactic (R a ) distances, the maximum 
(z max ) and minimum (z min ) distances to the Galactic plane, total orbital an- 
gular momentum (J z ), and the planar (e p ) and vertical (e^) orbital eccentrici- 
ties. Apogalactic and perigalactic distances are determined from the averaged 
maximum and minimum galactocentric distances of systems in the calculated 
Galactic orbits within the integration time of 3 Gyr, i.e. backwards in time 
over an interval of 3 Gyr. This integration time is chosen to correspond to 
12 or 15 revolutions around the Galactic center so that the averaged orbital 
parameters can be determined reliably. 



2. 3 Determination of Population Types 



In order to determine the population types of the CVs, we define some dy- 
namical and kinematical criteria. The criteria defined to select thick disc or 



5 



Table 1 

The data sample. The table can be obtained electronically. 



ID 


Star 


a (J2000.0) 


S (J2000.0) 




Ra 






Jz 




fin 






(hh:mm:ss) 


(»:':") 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc km s — 1 ) 






1 


AR And 


01:45:03.28 


+37:56:32.7 


7.894 


8.207 


-0.124 


0.124 


1790.1 


0.015 


0.019 


2 


LX And 


02:19:44.10 


+40:27:22.9 


6.517 


9.010 


-0.150 


0.150 


1681.0 


0.019 


0.160 


3 


PX And 


00:30:05.81 


+26:17:26.5 


7.001 


9.804 


-0.492 


0.491 


1815.7 


0.058 


0.167 


4 


RX And 


01:04:35.54 


+41:17:57.8 


7.984 


8.525 


-0.128 


0.128 


1837.4 


0.016 


0.033 


5 


V603 Aql 


18:48:54.64 


+00:35:02.9 


7.441 


7.957 


-0.081 


0.081 


1705.9 


0.010 


0.034 



halo CVs are as following: 



(1) 
(2) 

(3) 
(4) 



The total space velocity (Vug)- lAfgar et al.l ( 120121 ) showed that stars with 
V t ot > 100 km s _1 are thick disc or halo objects. 



The relative probab i lity fo r thick disc to thin disc membership (TD/D). 
Bensby et al.l (120031 . 120051 ) proposed that the stars with TD/D > 1 are 



thick disc or halo objects. 

The maximum distance from the Galactic plane (z max ). It is well known 
that the thick disk is the dominant component of the Galaxy between 1 
and 5 kpc above t he Galacti c disc while halo objects have z max values 
larger than 5 kpc ( Bilir et all 2008 ) . 

Vertical orbital eccentricity (e,,). iBilir et al. ( 2012 ) concluded that verti- 



cal eccentricities of Galactic orbits calculated for the thick-disc and halo 
stars are larger than ~0.1. 
(5) Planar eccentricit y of the Galactic orbit (e p ). Using main-sequence stars, 



Pauli et al.l (120031 ) found that stars with e p > 0.3 belong to the thick disc. 



TD/D is calculated using a purely kinematical approach desc ribed bylBensbv et al 



(120031 ). These probabilities were esti mated and lis t ed by lAk et al. 



According to the criteria defined by iBensby et al.l (120031 ) the stars are se- 
lected from four different TD/D intervals: TD/D < 0.1 (i.e. "high probabil- 
ity thin-disc stars"); 0.1 < TD/D < 1 (i.e. "low probability thin-disc stars"); 
1 < TD/D < 10 (i.e. "low probability thick-disc stars") and TD/D > 10 (i.e. 
"high probability thick-disc stars"). So, in this study a possible thick disc or 
halo CV is expected to fullfill at least two of these criteria. 



Using the criteria defined above, we find that nine of 159 CVs in the sample 
belong to thick disc population. The rest of the sample consists of thin disc 
systems. It is concluded that there are no halo CVs, as the maximum dis- 
tances to the Galactic plane z max of the CVs in the sample do not exceed 5 
kpc. Thick disc CVs found in this study are listed in Table 2 with their total 
space velocities (Vtot), with respect to LSR, types, orbital periods (P or b), arid 
the maximum distances to the Galactic plane (z max ). The relative probability 
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Table 2 

Probable thick disc CVs in the sample. Remarks indicate the fulfilled criteria as 
defi ned in the text. 

Criteria 









(1) 


(2) 


(3) 


(4) 


(5) 








Star 


Type 


Porb 


Vtot 


TD/D 


Z-m a x 


e„ 


e p 


Remark 






(hr) 


(km s — 1 ) 




(kpc) 












J2050-0536 


NL, AM 


2.299 


100.55 


14.5 


0.46 


0.038 


0.381 


1, 


2, 


5 


VV Pup 


NL, AM 


1.674 


151.63 


9602 


0.54 


0.036 


0.511 


1, 


2, 


5 


BY Cam 


NL, AM 


3.354 


142.74 


1064 


0.80 


0.060 


0.456 


1. 


2. 


5 


CT Ser 


N 


4.680 


80.58 


0.2 


1.00 


0.145 


0.243 


3, 


i 




V1043 Cen 


NL, AM 


4.190 


80.59 


8.3 


1.03 


0.102 


0.229 


2. 


3, 


4 


J1629+2635 


NL 


2.033 


72.11 


2.5 


1.04 


0.112 


0.184 


2. 


3, 


4 


IY UMa 


DN, SU 


1.774 


114.73 


70 


1.39 


0.148 


0.284 


1. 


2, 


3, 4 


J0813+4528 


DN, UG 


6.936 


78.34 


0.5 


1.77 


0.189 


0.198 


3, 


1 




AK Cnc 


DN, SU 


1.562 


154.72 


36015 


2.26 


0.195 


0.413 


1, 


2. 


3, 4, 5 



DN: dwarf novae, NL: nova-like stars, N: novae, SU: SU UMa type dwarf novae, UG: U Gem type dwarf novae, AM: AM 
Her (polars) type nova-like stars. 

0.6 



0.5 
0.4 

of 0.3 
0.2 
0.1 


50 100 150 200 

V tot (kms-i) 

Fig. 1. Vtot-e p diagram for the CVs in the sample. Thick disc and thin disc objects 
are indicated with filled and empty circles, respectively. Dashed lines represent limit 
values of criteria defined in the text. 

for the thick disc to thin disc membership (TD/D) is also listed in the Table 
2. e p and e v are planar and vertical eccentricities of the Galactic orbits, re- 
spectively. In Figs. 1-3, nine thick disc CVs are indicated in V tot -e p , e p -e v and 
e p -J z diagrams of the CV sample in this study, respectively. 

AK Cnc is the most probable thick disc CV in the sample. This object fulfills 
all the criteria we defined above. IY UMa is also one of the strongest thick 
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0.3 



0.2 



0.1 



a 



%g^Ma,o?^ 0m T ,„o, ,° i "I q o i 



0.1 0.2 0.3 0.4 0.5 O.I 



Fig. 2. e p -e v diagram for the CVs in the sample. Symbols and dashed lines are as 
in Fie. 1 
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o 2000 
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1500 
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0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 3. e p -J z diagram for the CVs in the sample. Symbols are as in Fig. 1 



disc candidates of CVs. IY UMa accomplishes four criteria but one: the planar 
eccentricity of its Galactic orbit is 0.284 which is very close to 0.3 (see, Fig. 
3).VV Pup, BY Cam and J2050-0536 accomplish three criteria: V tot > 100 km 
s -1 , TD/D > 1 and e p > 0.3 for these CVs (see, Fig. 2). For J1629+2635 and 



V1043 Cen, we found TD/D > 1, 



> 1 kpc and e v > 0.1. CT Ser and 



J0813+4528 match with two criteria. We found z max > 1 kpc and e v > 0.1 
for these two systems. AK Cnc, IY UMa, VV Pup, BY Cam a nd J2 050-0536 
were already suggested to be thick disc objects by lAk et al.l (120101 ) using a 
pure kinematical approach. Our results is consistent with this suggestion. 
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Fig. 4. Galactic orbits of the thick disc CVs found in this study projected onto 
X — Y and X — Z planes. Galactic orbits are calculated for an integration time of 
3 Gyr. 



3 Conclusion and discussions 



In this study we calculated, for the first time, the Galactic orbital parameters 
of a large number of CVs. We also determined their popula tion types using 



these calculations and the kinematical parameters as found by lAk et al.l (120101 ) 



Our analysis shows that almost all of the CVs in the sample are in the thin disc 
component of the Galaxy. Only nine CVs in the sample are thick disc stars. 
These objects are listed in Table 2. These are the CVs which cross the Galactic 
disc through their Galactic orbits. The Galactic orbits of these systems as 
projected on to X — Y and X — Z planes are shown in Fig. 4. Galactic orbits 
of the thick disc CVs in Fig. 4 were calculated for an integration time of 3 
Gyr, corresponding to 12-15 revolutions around the Galactic center. As the 
maximum distances from the Galactic plane (z max ) in Table 2 are not larger 
than 5 kpc, we conclude that there are no halo objects in the sample. 
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It is found that the fraction of thick disc CVs to thin disc CVs in the sample is 
~6% which is in very good agreement with th e fraction of Pop II field stars to 



Pop I field stars i n the Solar neighbourhood (IRobin et all Il996l ; iBuser et al. 



19991 ; iBilir et al.l . 120061 ) . This suggests that our CV sample is complete for 



the Solar neighbourhood. So, we conclude that statistical studies using this 
sample give reliable and self-consistent results. 



Ak et al.l (120101 ) found that the middle point of the period gap is at 2.62 hr 
for the objects in the sample. An inspection of the thick disc CVs in Table 
2 reveals that ~60% of thick disc CVs are located below the o r bital period 
gap. According to the population synthesis model of IStehle et al.l ( 119971 ). most 
Pop II CVs are expected to be found below the period gap. We conclude that 
the calculated fraction of thick disc CVs below the period gap to the thick 
disc systems above the period gap in this study is roughly consistent with 
the theoretically expected fraction. It is also interesting to note that ~45% of 
thick disc CVs found in this study are magnetic systems (AM Her stars). 

Some CVs have already been suggested to be the members of old populations 

i n the Galaxy. Object s prop o sed to be Pop II CV s inlHowell. Mitchell &: Warnock III! 

(ll987h . lDrissen et al.l (fliil . ISheets et al.l (120071) . lUthas et al.lfcoilh and llmamura fc Tanabe 
( 20121 ) are not included in our sample. Howell fc Szkody ( 1990 ) list 84 known 
or good candidates for being halo CVs by selecting high Galactic latitude ob- 
jects. Although 21 of them are included in our sample, only one of them (AK 
Cnc) is identified as a thick disc CV in this study. It must be emphasized 
that the Galactic latitude is not a reliable indicator for the population type 
of an object. Our analysis shows that the most probable thick disc CVs in 
Table 2 have high Galactic angular momentums, implying that they are in a 
population different than the thin disc (Fig. 3). As thick disc and halo objects 
have eccentric orbits, they can pass through the Solar neighbourhood during 
their nuclear evolution. That is why the analysis of Galactic orbits is very 
important in the determination of population types. 



The dispersion of the total space velocities is an indicator of population types, 
as the velocity dispersion of a group of objects is related to their kinematical 
age. The total s pace velocities o f the six systems with z max > 1 kpc in Table 2 
were taken from Ak et al. (l2010h and the dispersion of the t otal spa c e velo cities 



is calculated as 101 km s 1 . Using the equation given by IWielenl (119771 ) and 



Wielen et al.l (119921 ). the kinematical age of the six most probable thick disc 
CVs in the sample is calculated as 13 Gyr. Th is value is consistent with the 
age of the thick disc population in the Galaxy (IFeltzing fc Bensbyl . 120081 ) . 



Our study shows that there are thick disc CVs in the Solar neighbourhood. 
Our study also implies that halo CVs, if they are present, must be very rare 
in the vicinity of the Sun. Kinematical studies can reliably prove their pres- 
ence. However, current 7 velocity measurements are mostly dubious. That is 
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why we emphasize the importance of radial velocity studies of CVs. Further 
observational data can help to find more thick disc or halo CVs. Specially, 
investigations based on the data obtained from deep sky surveys could help to 
find more Pop II CVs. Detailed observations of these systems could give clues 
for their evolution. 
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